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(g) Procede d'enregistrement et de traitement d'une image et instrument de prise de vues pour la mise 
en oeuvre de ce proc£de. 

(5?) Procede d'enregistrement et de traitement 
d'une image d'une scene situee dans un plan P 
donne rapporte a deux directions perpendlcu- 
laires x, y, tadite image etant foumie par un 
instrument (100) de prise de vues presentant, 
dans la direction u (u = x, y) du plan P, une 
fonction de transfert de modulation H u (f) de 
frequence de coupure f^, f designant la fre- 
quence spatiale, procede comprenant, d'une 
part, une premiere etape d'enregistrement de 
ladite image a une frequence d'echantillonnage 
f eu de I'image dans la direction u, egale a kf cu , k 
etant un nombre entier superieur a 1 et, d'autre 
part, une deuxieme etape de traitement de I'i- 
mage enregistree consistant a appliquer a ladite 
image un flltre de deconvolution. 

Selon Tinvention, ledit flltre de deconvolution 
est egal, dans le domaine frequentiel et pour 
chaque direction u, au produit de I'inverse de la 
fonction de transfert de modulation H u (f) par 
une fonction de troncature l u (f) ajustable. 

Application a I'observation de la Terre a partir 
de satellites artificiels. 
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La pr6sente invention concerne un procede d'enregistrement et de traitement d'une image d'une scene 
situ6e dans un plan P donne rapporte a deux directions perpendiculaires x, y, ladite image etant fournie par 
un instrument de prise de vues pr6sentant, dans la direction u (u = x. y) du plan P, une fonction de transfert 
de modulation Hu (f) de frequence de coupure f cu , f designant la frequence spatiale. 
5 L'invention trouve une application particulierement avantageuse dans le domaine de I'observation de la 

Terre a partirde satellites artificiels. 

Les instruments de prise de vues utilises dans les systemes satellitaires d'observation de la Terre compor- 
tent, soit un detecteur eiementaire dont I'image dans le plan P parcourt la scene a enregistrer a I'aide d'un dis- 
positif de balayage m6canique couple au mouvementde translation ou de rotation du satellite, soit une barrette 
10 Iin6aire de detecteurs eiementaires, par exemple du type a couplage de charges (CCD), avec balayage en 
translation dans la direction perpendiculaire a I'axe de ladite barrette (balayage dit "push broom"). 

De fac on usuelle, le proc6d6 conforme au preambule comprend une premiere etape d'enregistrement de 
rimage a une frequence d'6chantillonnage f ou dans la direction u. egale a la frequence de coupure f cu dont Tin- 
verse correspond a la resolution spatiale au niveau du plan P d observation. Cette resolution est en g6n6ral 
is Iimit6e par les dimensions du (des) detecteur(s) eiementaire(s) utilises). 

Toutefois, du fait que le dispositif de formation d'images de imstrument de prise de vues, un telescope 
par exemple, n'est pas lui-meme parfait, la fonction de transfert de modulation H u (f) n'est pas egale a 1 sur 
tout le spectre de frequence, ce qui conduit a une degradation de I'image enregistree de la scene a observer. 
II est possible de traiter cette image en iui appliquant un f iltre de deconvolution. Si on appeile P u (f) et S u (f) les 
20 spectres de frequence respectifs dans la direction u de la scene a observer et de I'image enregistree, on a la 
relation : 

SJ!) = H u (f) P u (0 

Pour restituer P u (f), il faut done appliquer a S u (f) un f iltre de deconvolution qui dans le domaine fr6quentiel est 
6gal a I'inverse de la fonction de transfert de modulation 

Pu(f) = Hio Su(0 

Ce qui s'exprrme dans le domaine spatial par le produit de convolution : 

P(u) = TF " 1 t^)l® s ( u > 

30 TF 1 d6signant la transformee de Fourier inverse, et p(u) et s(u) etant les fonctions spatiales repr6sentant res- 
pectivement la sc6ne et son image. 

Dependant la th6orie du traitement du signal, notamment le theoreme de Shannon, montre que seules les 
frequences inferieures ou 6gales a la moitie de la frequence d'echantillonnage peuvent etre restitu6es par de- 
convolution, cette limitation 6tant due au ph6nomene general de repliement spectral. De plus, lors de la de- 

35 convolution, les hautes frequences, comprises entre f ou /2 et f eu , viennent perturber les basses frequences 
comprises entre 0 et f eu / 2 - 

II a 6te propose dans Tarticfe "Repliement spectral lie a I'utilisation de detecteurs constitues de cellules 
CCD" de G. Jacquemod etal., paru dans la revue 'Traitement du Signal" Vol. 5, n° 1, 1988, de limiterl'influence 
du repliement spectral par un sur6chantillonnage de I'image, la frequence d'echantillonnage f eu etant prise 6ga- 
^0 le a kf cu , k etant un nombre entier sup6rieur a 1. Dans ce cas, il est possible de restituer les frequences au 
moins jusqu'a la frequence de coupure f^. D'autre part, comme au-dela de f cu fa fonction de transfert de mo- 
dulation est pratiquement nulle, les frequences superieures a f cu sont sans effet sur les basses frequences in- 
ferieures a f cu . 

On dispose done d'un moyen d'augmenter d'un facteur k la resolution dans la direction u d'un instrument 
45 de prise de vues donn6 ou, selon un autre aspect du procede, de r6aliser un instrument de prise de vues ayant 
les m£mes performances de resolution mais de moindre encombrement dans la mesure ou la distance focale 
du dispositif de formation d'images peut etre divis6e par k tout en conservant la meme ouverture numerique, 
et done le m§me rapport signal a bruit. Cette derniere possibilite pr6sente un avantage tr6s important dans le 
cas de I'application spatiale envisag6e ci-dessus puisque ne serait-ce qu'un facteur k = 2 permettrait un gain 
so de 8 en volume et 4 a 5 sur la masse de Instrument de prise de vues. d'ou une diminution appreciable de la 
charge utile du satellite. 

II faut cependant signaler que les avantages pr6cites ne sont pas sans inconv6nients. En effet, que ce 
soit dans le cas d'une augmentation de la resolution d'un instrument de prise de vues donn6, ou dans celui 
de I'utilisation d'un dispositif de formation d'images de distance focale r6duite, I'op6ration de deconvolution a 
55 pour consequence une diminution du rapport signal a bruit, ceci etant du au fait que, compte tenu du spectre 
de bruit B u (f), le signal radiometrique S u (f) de I'image s'ecrit : 

S u (0 = H u (f) P u (f) + B u (f). 
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[.'application du f iltre de deconvolution conduit & une scene restitute : 

u( ° H u (f) Pu(f) + H U (D 
La fonction de transfert de modulation H u (f) etant au plus egale & 1 , on comprend que I'effet du bruit B u (f) 
5 a tendance & augmenter. Aussi, le probleme technique & resoudre par I'objet de la presente invention est de 
proposer un procede conforme au preambule, comprenant, d'une part, une premiere etape d ? enregistrement 
de ladite image d une frequence d'echantillonnage f eu de I'image dans la direction u egale & kf CUf k etant un 
nombre entier superieur & 1 , et, d'autre part, une deuxieme etape de traitement de l'image enregistree consis- 
tant £ appliquer & ladite image un f iltre de deconvolution, procede qui permettrait, lors de I'operation de de- 
convolution, de conserver I'essentiel des avantages du procede de surechantillonnage tout en ayant la possi- 
bility de toujours maitriser le niveau de bruit final a une valeur acceptable. 

La solution au probleme technique pose consiste, selon la presente invention, en ce que ledit f iltre de de- 
convolution est egal, dans le domaine frequentiel et pour chaque direction u, au produit de I'inverse de la fonc- 
tion de transfert de modulation H u (f) par une fonction de troncature l u (f) ajustable. 
*5 Ainsi, la scene reconstitute a pour spectre de fr6quence : 

Pu'W = l„(0 P u (0 + l u (f) ~^ 

On notera que l effet de la fonction de troncature l u (f) varie en sens contraire selon que Ton considere la 
qualite image ou le rapport signal d bruit. Si Ton desire limiter sensiblement le bruit final, on choisira une fonc- 
tion de troncature faible, au detriment de la resolution spatiale representee par l u (f) P u (f), ce qui n'est pas tres 
penalisant dans le cas d'un instrument de prise de vues donne dontle surechantillonnage augmente la reso- 
lution et reduit le signal. Par contre, pour un instrument sous-dimensionne dont le surechantillonnage permet 
de maintenir la resolution avec une meme ouverture et done le meme signal, on pourra choisir une fonction 
de troncature l u (f) plus importante. Bien entendu, selon 1'instrument considere, la fonction de troncature sera 
choisie et adaptee, en tenant compte des caracteristiques de la mission & accomplir. 

La description qui va suivre en regard des dessins annexes, donnes £ titre d'exemples non limitatifs fera 
bien comprendre en quoi consiste invention et comment elle peut etre realisee. 

La figure 1 est un schema d un instrument de prise de vues pour la mise en oeuvre du procede selon In- 
vention. 

La figure 2 montre une fonction de transfert de modulation typique de I'instrument de la figure 1. 
La figure 3 montre une fonction de troncature applicable en combinaison avec la fonction de transfert de 
modulation de la figure 2. 

La figure 4a montre un mode de realisation du surechantillonnage colonne pour une barrette Hneaire de 
detecteurs. 

La figure 4b montre un mode de realisation du surechantillonnage ligne pour une barrette lineaire de de- 
tecteurs. 

La figure 5 est un premier montage de barrettes lineaires realisant un surechantillonnage ligne d'ordre 2. 
La figure 6 est un deuxieme montage de barrettes lineaires realisant un surechantillonnage ligne d'ordre 

2. 

La figure 1 represente de facon schematique un instrument 100 de prise de vues, installe par exemple £ 
bord d*un satellite artif iciel de la Terre, et destine d fournir une image d'une scene situee dans un plan P donne, 
rapporte £ deux directions perpendiculaires x, y. Dans rapplication envisagee, le plan P correspond au sol. 

L'instrument 100 de prise de vues est constitu6, d'une part, d'un dispositif de formation d'images sche- 
matise par une lentille 11 0, mais qui pourrait etre tout aussi bien un telescope & miroirs par exemple, et, d'autre 
part, d'une barrette lineaire 120 de detecteurs elementaires 121 du type CCD, placee au foyer F de la lentille 
110. Dans la suite de la description et sans que cela ne soit limitatif, les detecteurs elementaires 121 sont sup- 
poses de forme carree. Comme on peut le voir sur la figure 1 , le mode de balayage du plan P est du type "pus- 
hbroom" dans lequel I'image d'axe x de la barrette lineaire 120 se deplace dans la direction perpendiculaire y 
avec une vitesse V egale £ la vitesse au sol du satellite. 

L'instrument 100 de prise de vues presente dans chaque direction u = x, y du plan P une fonction de trans- 
fert de modulation. H u (f) dont une allure typique est donnee sur la figure 2. Cette fonction de transfert de mo- 
dulation s'annule pour une frequence de coupure f^ correspondant a la resolution de l'instrument 100. Gene- 
ralement, le pouvoir de resolution est limite par la structure discontinue des detecteurs plus que par la qualite 
du dispositif de formation d'images. Si dx represente la taille des detecteurs elementaires 121, la resolution 
au sol Au dans la direction u du plan P est donnee par 

Au = tLdx 
fo 
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H etant I'altitude du satellite et fo la distance focale du dispositif de formation d'images. A titre d'exemple, pour 
un satellite dont I'orbite est situ6e a 800 km de la Terre, une distance focale de 1 m et des detecteurs eie'men- 
taires de 7 u.m de c6te, on atteint une resolution au sol de 5,6 m. 
La frequence de coupure f^ est alors sensiblement 6gale a 

Tcu Au*Hdx 

L'image de la sc6ne a observer est d'abord enregistr6e a une frequence d'6chantillonnage f eu dans la di- 
rection u egale a kf^, k etant un nombre entier superieur a 1. Dans le mode de realisation particulier montre 
a la figure 2, le coefficient de sur6chantilIonnage k a ete pris 6gal a 2, de sorte que f ou vaut 

f = Of = 

•eu ^ 'cu . 

AU , 

Aprfcs cette premiere etape d'enregistrement, rimage enregistr6e S u (f) est soumise a une operation de trai- 
tement consistant a lut appliquer un f Htre de d6convolution 6gal, dans le domaine frequentiel et dans chaque 
direction u t au produit de I'inverse de la fonction de transfert de modulation H u (f) par une fonction de troncature 

15 lu(f). 

La figure 3 montre un exemple d'une telle fonction de troncature analytiquement donnee par 



1 si f $i { 



I u (f) J cos 2 71 (f - f l } si t l4 t <f cu 



2(f cu- f l> 



•0 si f * f ( 



Cette courbe est connue sous le nom de fenStre de Hanning. 

Comme cela a et6 expliqu6 plus haut, la fonction de troncature l u (f) est ajustable de manidre a obtenir, en 
30 fonction de la mission envisag6e, le meilleur compromis entre resolution spatiale et rapport signal a bruit. 

Les figures 4a et 4b indiquent respectivement comment, dans le cas d'une barrette Iin6aire de CCD, un 
sur6chantillonnage par 2 en colonne et en ligne peut dtre obtenu. 

En ce qui concerne le surechantillonnage colonne (Figure 4a), it suff it de disposer d'un temps d'integration 
T, de la barrette 6gal au temps At qui met la trace au sol du satellite pour parcourir la distance A y/2 soit : 

Ti = A t = = — 
2V 2V 

Pour realiser le surechantillonnage ligne (Figure 4b), ii suff it de disposer de deux barrettes eiementaires 
paralieies. de mfime pas dx, d6caiees Tune par rapport a I'autre de la distance dx/2 de sorte que les images 
de ces deux barrettes dans le plan P soient d£cal£es de A x/2 = A u/2. 

40 La figure 5 donne un premier montage de deux barrettes lin£aires realisant uri surechantillonnage ligne 

d'ordre 2. Les barrettes 120* et 120" sont montees chacune sur une face d'un prisme respectif 210', 210", les- 
dits prismes etant acco!6s de facon a constituer un diviseur optique 200 de ligne. Le flux incident 300 se partage 
a ['interface 220 des deux prismes formant lame semi-r6fiechissante. Chaque barrette Iin6aire 120*et 120" ne 
recoit done que la moitie de l'6nergie lumineuse provenant de la scene a observer, d'ou une limitation du rapport 

45 signal a bruit. 

Le montage represents sur la figure 6 permet d'obtenir un sur6chantillonnage d'ordre 2 avec un signal mul- 
tiplie par un facteur 2 par rapport au montage precedent. Les deux barrettes Iin6aires 120', 120", d6cal6es de 
dx/2, sont r6alis6es dans un mSme composant 400, et, de facon a assurer la synchronisation des images en- 
registr6es s6par6ment par chacune des barrettes, le retard a I'enregistrement de I'une des barrettes parrap- 
50 port a Tautre doit etre egat a k'T,, k' etant un nombre entier non nul, ce qui impose que la distance d s6parant 
les barrettes soit 6gale a k'dx/2. Dans ce deuxi6me montage, chaque ligne de detecteurs recoit I'int6gralite du 
flux lumineux incident, d'oii un signal double de celui correspondant au montage de la figure 5. Dans le cas 
ou le bruit photonique est preponderant, le rapport signal a bruit est multiplie par un facteur 6gal a V2. 

II faut egalement noter que le dernier montage peut se g6n6raliser au cas d'un sur6chantillonnage d'ordre 
55 k superieur a 2. II suff it de disposer dans un m£me composant de k barrettes Iin6aires decaiees de dx/k et 
distantes de k'dx/k. Le nombre k est seulement limite par les possibilites technologiques de r6aliser un compo- 
sant a plusieurs barrettes. 

Enf in. I'invention a 6t6 d6crite en mettant en oeuvre une plurality de detecteurs eiementaires assembles 
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en barrettes, mais i! est bien entendu qu'elle pourrait tout aussi bien §tre r6alis6e avec un seul d6tecteur 6I6- 
mentaire. Dans ce cas, I'enregistrement de I'image sur le d6tecteur est effectu6 par balayage de Timage dudit 
d6tecteur 6l6mentaire dans les deux directions du plan P avec sur6chantillonnage sur chacun des balayages. 



Revendications 

1. Proc6d6 d'enregistrement et de traitement d'une image d'une sc6ne situ6e dans un plan P donn6 rapport6 
a deux directions perpendiculaires x, y, ladite image 6tant fournie par un instrument (1 00) de prise de vues 
pr6sentant, dans la direction u (u = x, y) du plan P, une fonction de transfert de modulation H u (f) de fr6- 
quence de coupure f^, f d6signant la frequence spatiale, proc6d6 comprenant, d'une part, une premiere 
6tape d'enregistrement de ladite image & une fr6quence d*6chantillonnage de I'image dans la direction 
u, 6gale £ kf cu , k 6tant un nombre entier supSrieur & 1 et, d'autre part, une deuxteme 6tape de traitement 
de I'image enregistr6e consistant d appliquer & ladite image un f iltre de d6convolution, caract6ris6 en ce 
que ledit f iltre de d6convolution est 6gal, dans le domaine fr6quentiel et pour chaque direction u, au produit 
de rinverse de la fonction de transfert de modulation H u (f) par une fonction de troncature l u (f) ajustable. 

2. Proc6d6 selon la revendicatipn 1, caract6ris6 en ce que, ledit instrument de prise de vues comportant au 
moins un d6tecteur 6l6mentaire (121) et la fr6quence de coupure f cu 6tant sensiblement 6gale & I'inverse 
de la dimension Au de I'image dudit d6tecteur 6l6mentaire dans la direction u du plan P, la frequence 
d'6chantillonnage f eu est 6gale d k/ Au. 

3. Proc6d6 selon I'une des revendications 1 ou 2, caracteris6 en ce que, ledit instrument de prise de vues 
comportant un d6tecteur 6l6mentaire (121), I'enregistrement de ladite image est r6alis6 par balayage de 
I'image dudit d6tecteur 6!6mentaire dans les deux directions x, y du plan P. 

4. Proc6d6 selon I'une des revendications 1 ou 2, caract6ris6 en ce que, ledit instrument (100) de prise de 
vues comportant au moins une barrette Iin6aire (120) d'une plurality de detecteurs 6!6mentaires (121) 
dont les images dans le plan P sont align6es dans la direction x, Penregistrement de ladite image est r6alis6 
par balayage de I'image de ladite barrette Iin6aire (120) dans la direction y. 

5. Instrument de prise de vues pour la mise en oeuvre du proc6d6 selon la revendication 4, caract6ris6 en 
ce que ledit instrument (100) comporte k barrettes Iin6aires paraJI6!es (120, 120"), de mSme pas dx, d6- 
cal6es successivement de la distance dx/k et s6par6es I'une de I'autre de la distance k'dx/k, k' 6tant un 
nombre entier non nul. 

6. Instrument de prise de vues pour la mise en oeuvre du proc6d6 selon la revendication 4, caract6ris6 en 
ce que, le nombre k 6tant 6gal & 2, ledit instrument (100) comporte au moins deux barrettes Iin6aires pa- 
rallfeles (120', 120") de m£me pas dx, disposes sur un diviseur optique de ligne et d6cal6es I'une par 
rapport a I'autre de la distance dx/2. 
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FIG. 4a 
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FIG.4b 
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